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ABSTRACT 



The cosmic star formation histories are evaluated for different minimum masses of 
the initial halo structures, with allowance for realistic gas outflows. With a minimum 
halo mass of 10 7 -10 8 Mq and a moderate outflow efficiency, we reproduce both the 
current baryon fraction and the early chemical enrichment of the IGM. The intensity 
of the formation rate of "normal" stars is also well constrained by the observations: it 
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has to be dominated by star formation in elliptical galaxies, except perhaps at very low 
redshift. The fraction of baryons in stars is predicted as are also the type la and II 
supernova event rates. Comparison with SN observations in the redshift range z = — 2 
allows us to set strong constraints on the time delay of type la supernovae (a total delay 
of ~4 Gyr is required to fit the data) , the lower end of the mass range of the progenitors 
(2-8 M Q ) and the fraction of white dwarfs that reproduce the type la supernova (about 
1 per cent). The intensity in the initial starburst of zero metallicity stars below 270 M Q 
must be limited in order to avoid premature overenrichment of the IGM. Only about 
10 - 20 % of the metals present in the IGM at z = have been produced by population 
III stars at very high z. The remaining 80 - 90 % are ejected later by galaxies forming 
normal stars, with a maximum outflow efficiency occurring at a redshift of about 5. We 
conclude that 10~ 3 of the mass in baryons must lie in first massive stars in order to 
produce enough ionizing photons to allow early reionization of the IGM by z ~ 15. 

Subject headings: Cosmology: theory — nuclear reactions, nucleosynthesis, abundances 
— stars: evolution — Galaxy: evolution 

1. Introduction 

Chemical evolution is a key to understanding the hybrid role of massive star formation in the 
early universe with regard to the epoch of reionization, the heavy element abundances of the oldest 
stars and the high z intergalactic medium, and the mass outflows associated with galaxy formation. 
In addition, predicted supernova rates provide us with an independent probe of the early epoch 
of star formation. Combining abundance and supernova rate predictions allows us to develop an 
improved understanding of both the cosmic star formation history and of the enrichment of the 
IGM, as well as to elucidate the nature of Population III (see e.g. the recent review of Ciardi & 
Ferrara (2005)). 

The prevalent view is that the first stars encompassed the mass range 100 to 1000 M . Recent 
support for this possibility stems from the need to reionize the Universe at high redshift (Cen 2003; 
Haiman & Holder 2003; Wyithe & Loeb 2003b; Bromm 2004) as indicated by the WMAP first-year 
data (Kogut et al. 2003). The possibility for early reionization by the first galaxies was considered 
by Ciardi et al. (2003) and further support for this hypothesis is based on chemical abundance 
patterns at low metallicity (Wasserburg & Qian 2000; Oh et al. 2001; Qian & Wasserburg 2001, 
2005). 

However, the robustness of the conclusion that very massive stars (VMS) were necessarily 
present among the first stars has been questioned (Venkatesan & Truran 2003; Tumlinson et al. 
2004; Tumlinson 2005). A "normal" initial mass function (IMF) may be capable of producing 
ionization consistent with WMAP (Venkatesan et al. 2003; Wyithe & Loeb 2003b), and within 
these constraints, it was argued (Venkatesan et al. 2004) that a broad set of chemical abundances 
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may be better fit using the yields of Umeda & Nomoto (2003) for stars with masses in the 1- 50 
M range than with the yields from pair-instabilty supernovae (PISN) (Heger & Woosley 2002). 
Indeed, these results were confirmed in Daigne et al. (2004, hereafter DOSVA) where it was argued 
that a top-heavy IMF without VMS supplied a better fit to low metallicity abundance data while 
still accounting for the early re-ionization of the Universe when using a detailed model of cosmic 
chemical evolution. Furthermore, when VMS were assumed to be the predominant source for Pop 
III, it was found that there could not be enough of these stars to have accounted for reionization 
without producing anomalous chemical abundance ratios in both old halo stars and in the IGM. 
Even partial suppression of the yields of pair instability SNe via rotation and fall-back should not 
vitiate this conclusion in view of the anomalous abundance ratios and the large discrepancy found 
by DOSVA. 

In this paper, we improve on an earlier model for primordial star formation and chemical 
evolution (DOVSA) by using a more realistic model of structure formation to study the role of 
baryonic outflows and the enrichment of the IGM. Our principal result is that only if Population 
III predominantly consisted of stars in the "normal" mass range, but with a top-heavy IMF, that 
is to say dominated by stars in the mass range 40 to 100 M , can one simultaneously account for 
reionization and chemical enrichment, both of galaxies and of the IGM. However, this mass range 
is not capable of sufficient Si production and therefore, more massive stars may be necessary to 
account for the inferred abundance of Si in the IGM. Such a bimodal approach to star formation 
allows us to account both for the cosmic star formation history and global chemical evolution of 
the universe. The predicted supernova rates agree with those observed. We also infer a baryonic 
outflow rate that allows us to reconcile the observed baryon fraction in massive galaxies with that 
observed in clusters, the CMB and the Lyman alpha forest, as well as that predicted by primordial 
nucleosynthesis . 

The outline of this paper is as follows. In § 2 we present the details of the chemical evolution 
model, especially the treatment of the structure formation and of the galactic winds. We also 
outline the key parameters used in the models which affect our conclusions. In § 3, we discuss the 
normal mode of star formation and we show that the minimum masses of the halos of star forming 
structures as well as the efficiency of the galactic winds are well constrained by the observations. 
This allows us to select a few realistic models. In § 4, we compute the expected SNII and SNIa 
rates in these models and deduce the time delay inferred for SNIa to fit the data. § 5 is devoted to 
the addition of an initial starburst of massive stars in our scenario. We summarize our results in 
5 6. 
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2. Hierarchical formation scenario and associated galactic winds 

2.1. A cosmic evolutionary model 

The chemical evolution model used in this paper has been described in DOVSA. It is a gener- 
alized version of standard models designed to follow one specific structure such as the Milky Way 
(for a review, see Tinsley (1980)). We describe baryons in the Universe by two large reservoirs. The 
first is associated with collapsed structures (hereafter the "structures") and is divided in two sub- 
reservoirs: the gas (hereafter the "interstellar medium" or ISM) and the stars and their remnants 
(hereafter the "stars"). The second reservoir corresponds to the medium in between the collapsed 
structures (hereafter the "intergalactic medium" or IGM). The evolution of the baryonic mass of 
these reservoirs, i.e. Migm(*) of the IGM, Mism(*) of the ISM and M*(i) of the stars is governed 
by a set of differential equations (see section 2 in DOVSA): 

dMi G M dM stTUCt 

= _ = -a b (f) + (f), (1) 



dM * , T ,m u\ a dM ™ M dM strnct dM* 

— : — = Wit) — e(t) and = — . 2) 

dt w w dt dt dt v ' 

In addition, we have Mism(^) + M*(t) = M struct (t), corresponding to the total baryonic mass of the 
structures, and Migm(^) + M stTUCi (t) = constant, which is the total baryonic mass of the Universe. 
As can be seen, these equations are controlled by four rates which represent four fundamental 
processes (see sketch in Figure 1): the formation of structures through the accretion of baryons 
from the IGM, ab(i); the formation of stars through the transfer of baryons from the ISM, w(t); 
the ejection of enriched gas by stars, e(t) and the outflow of baryons from the structures into the 
IGM, o(t). 

We track the chemical composition of the ISM and the IGM separately as a function of time (or 
redshift). The differential equations governing the evolution of the mass fraction Xj SM (Xj GM ) of 
element i in the ISM (IGM) are given by equations (6) and (7) in section 2 in DOVSA. In addition, 
we also compute the ionizing UV fluxes from the stars (equation (12) in section 3 in DOVSA) and 
the rate of explosive events (type la and gravitational collapse supernovae). 

The age t of the Universe which appears in the equations is related to the redshift by 

dt _ 9.78k- 1 Gyr 

~d~z ~ / =z ' W 

(i + z)^n A + n m (i + z ) 3 

assuming the cosmological parameters of the so-called "concordance model", with a density of 
matter Q m = 0.27 and a density of "dark energy" = 0.73 (Spergel et al. 2003), and taking 
Hq = 71 km/s/Mpc for the Hubble constant (h = 0.71). This allows us to trace all the quantities we 
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Fig. 1. — A sketch of our cosmic evolutionary model. The three baryon reservoirs (IGM, ISM and 
stars) exchange mass via four physical processes (structure formation, star formation, gas 
ejection by stars, e, outflow from the structures, o). 
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describe as a function of redshift. The input stellar data (lifetimes, mass and type of the remnant, 
metal yields, UV flux) are taken to be dependent both of the mass and the metallicity of the star 
(see DOVSA for more details). 

Building upon our original model, we have included the following improvements: (i) the baryon 
accretion rate, a^(t), by the structures is now computed in the framework of the hierarchical scenario 
of structure formation, (Press & Schechter 1974; Sheth k, Tormen 1999; Jenkins et al. 2001; Wyithe 
& Loeb 2003a). This is described in section 2.2. (ii) the baryon outflow rate, o(t), from the 
structures now includes a redshift-dependent efficiency, which accounts for the increasing escape 
velocity of the structures as the galaxy assembly is in progress. This is described in section 2.3. 
Note that in DOVSA, the outflow contained two sources: o{t) = o w (i) + osn(^) where o w (i) is 
a global outflow powered by stellar explosions (galactic winds) and osn(^) corresponds to stellar 
supernova ejecta that are flushed directly out of the structures. Since it was shown in DOVSA 
that the second term has a very small effect on the results, we have neglected it in this work (in 
DOVSA's notation, this is equivalent to setting a = 0). 



As seen in DOVSA, the term a^^t) which stands for the process of structure formation can 
have a strong impact on the high redshift evolution, intensity of the star formation process, ionizing 
flux, and the enrichment of the IGM, as it governs the size of the reservoir of baryons available for 
star formation. In DOVSA, although this term was estimated in an oversimplified way, we assumed 
only that structure formation efficiency decays exponentially. This first attempt allowed us to point 
out the key role of this term in the generalized chemical evolution model we have implemented. 
Here, we improve the model with a more realistic description of structure formation. 

We adopt the framework of the hierarchical scenario where small structures are formed first. 
At redshift z, the comoving density of dark matter halos in the mass range [M, M + dM] is 
/ PS (M,z)dM, with 



where /9dm is the comoving dark matter density. The distribution function of halos fps (M, z) is 
computed using the method described in Jenkins et al. (2001) using a code provided by A. Jenkins. 
It follows the standard theory (Press &: Schechter 1974), including the modification of Sheth &: 
Tormen (1999) and assumes a primordial power spectrum with a power-law index n = 1 and the 
fitting formula to the exact transfer function for non-baryonic cold dark matter given by Bond & 
Efstathiou (1984). We adopt a rms amplitude a% = 0.9 for mass density fluctuations in a sphere of 
radius 8 hr x Mpc. 



2.2. 



Hierarchical formation scenario 




(4) 



- 7- 



We assume that the baryon distribution traces the dark matter distribution without any bias 
so that the density of baryons is just proportional to the density of dark matter by a factor 
Qft/ (fi m — f2&). We take a baryonic density f2b = 0.044 (Spergel et al. 2003). We parametrize 
the fact that stars can form only in structures which are suitably dense by defining the minimum 
mass M min of a dark matter halo of the collapsed structures where star formation occurs. This mass 
could be related in principle with the critical temperature at which the cooling processes become 
efficient enough to allow star formation. In fact, this critical temperature, and hence the minimum 
mass M m in, should evolve with redshift z as the cooling processes of the hot gas in structures 
depend strongly on the chemical composition and ionizing state of the gas. It is however beyond 
the scope of this study to include such a detailed analysis so we prefer to keep M m \ n constant and 
to consider it as a free parameter of the model. The fraction of baryons at redshift z which are in 
such structures is then given by 



fb. 



struct 



(z) = 



f™dM Mf PS (M,z) 
Therefore, the mass flux can be estimated by 



(5) 



a h (t) 
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dz 



= 1.2h 3 M /yr/Mpc 3 (J^j (1 + z) yJn A + Sl m (l + . 



dfb, 



struct 



dz 



(6) 



This is the new expression of a\,(t) used in the model. Even if it represents a significative im- 
provement compared to the treatment of the same term in DOVSA, one should keep in mind that 
this derivation is only partially consistent, as the actual fraction of baryons in collapsed struc- 
tures with dark matter halos of mass M > M nun will differ from the value given by equation 5 
due to the outflow of baryons from the star-forming structures towards the IGM. However, since 
o(t) is always small compared to a^t), this is a small correction for most of the evolutionary history. 



Finally, this formalism allows us to indirectly fix the initial redshift z; n it where the first stars 
form in the simulation. This redshift is a priori smaller than the redshift zps(M nun ) at which the 
first dark matter halos of mass M m i n appear, as the star formation process can be delayed due to 
the time necessary to cool the collapsing gas. Instead of fixing the value of z- m i t , we prefer to fix 
the size of the gas reservoir when the star formation starts, i.e. the initial fraction /b, s truct(^imt) 
of baryons in collapsed structures which form stars. We usually adopt instruct (^init) = 1 % but we 
also consider 0.1 % and 5 % for some specific cases. Supernova feedback generically gives a star 
formation efficiency of 1 or 2 percent per dynamical time (Silk 2003) and an inefficiency of this 
order is inferred in essentially all current epoch star- forming galaxies (Kennicutt 1998). 

The corresponding initial redshift Zi n a is computed from equation 5 for a given minimum mass 
-^min (see table 1). Note that this redshift is usually smaller than that found in hydro dynamical 
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simulations of collapsing gas for the redshift of formation of the first zero-metallicity stars (popula- 
tion III stars): z ~ 20 — 30 (Bromm & Larson 2004; Yoshida et al. 2004; Bromm 2004). However, 
the first stars formed in these simulations are found to form in collapsed structures with dark 
matter halos of mass M ~ 10 6 M Q which represents only a very small fraction (< 10~ 3 ) of the 
baryons in the Universe at these redshifts. Therefore, this very primordial star formation activity 
is completely negligible in comparison with the global evolution of the Universe. 



2.3. Outflows 

We compute the outflow by 

2e f mup 

°(*) = ~1~T~\ / dm ®( m )^ (* - r ( m )) ^kin(m) , (7) 

v esc\ z ) ./max (8 M ,m d (t)) 

where 3>(m) is the initial mass function (IMF) of stars, r(m) is the lifetime of a star of mass m and 
is the mass of stars that die at age t, E^ n (m) is the kinetic energy released by the explosion of 
a star of mass m, e is the fraction of the kinetic energy of supernovae that is available to power the 
outflow and v^ sc (z) is the mean square of the escape velocity of structures at redshift z (see e.g. 
Scully et al. 1997). The escape velocity is obtained by mass-averaging the gravitational potential 
over the mass distribution at redshift z: 

!r im JMf^M,z)M(2GM/R) 
eSc[ > Jm^JM fi*{M,z)M ' W 

where R is the radius of a dark matter halo of mass M. The resulting velocity is plotted in Figure 2 
for different minimum masses M mm . Notice the strong evolution with z: outflows at high redshift 
are much more efficient and a large fraction (if not all) of the gas in the structures is ejected and 
mixed in the intergalactic medium, thus contributing to the early enrichment of the IGM. 

The total efficiency of this outflow process will depend both on e and v^ sc (z) and will therefore 
vary with redshift. An approximation of the efficiency can be obtained in the "instantaneous 
recycling limit" where the lifetime of massive stars is neglected by comparing the instantaneous 
outflow rate with the star formation rate. In this case we have 

_ o(f)_ ^ 2e e kin 

^Outflow — lT ,/j\ — O f \1 \") 

where ekin is the kinetic energy per stellar mass, averaged over the initial mass function: 

rm up 

ekin = / dm <$(m)E kin (m) . (10) 

Jm ini 

For a power-law IMF of slope — (1 + x), a minimal mass m in f and a constant supernova kinetic 
energy E^i n (m) ~ constant, this kinetic energy per stellar mass can be approximated by 

x — 1 .Ekin j (8 M Q /m in f)~ x if m in f < 8 M (normal mode) . . 

e^n — x s • (1-U 

x m m f | 1 if m-inf > 8 M (massive mode) 
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Fig. 2. — The averaged escape velocity from structures, v esc , as a function of redshift is shown for 
two choices of M m j n . The solid curves correspond to the velocity in units of km/s as indicated by 
the scale on the left, whereas, the dotted curves correspond to v 2 in erg/M with values indicated 
on the right. 
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The derivation of eq. 11 assumes a properly normalized IMF (see eq. 16 below). Typical values for 
v esc( z ) & t different redshifts can be obtained from Figure 2. These result in the following efficiencies 
for x = 1.3 and M min = 10 7 M : 

r — * 0052 (m) (ife) (oSfeF ( »xi££UJ ~' < 12> 

for a normal mode of stellar formation at z ~ 1, 



^outflow 



9 / f W gkin \ / m inf \ a3 / ^g sc (2) \ 
' V0.01 A 10 51 erg / V 0.1 M m / I 8 x 10 46 erg/M m / 1 ' 



l 5i erg 7 VO.l M y V 8 x 1046 erg/M Q 
for a normal mode of stellar formation at z ~ 8, 

#kin \ / m inf / Vesc(^) 



^outflow - 5.8 ( Q0l ) (^ 1051 k " rg j (^ 4() ^ q J ly 2 x 1O 46 C erg/M0 J ( 14 ) 

for a massive mode of stellar formation at z ~ 16 (m in f = 40 M ) and 

routfiow - 8-2 (— ) ( 1Q k 5i ) ( 140 M m ) (2 x 1046 C( eri/M m ) (15) 



0.01/ V5 x 10 51 ergy Vl40 M y V 2 x !0 46 erg/M 

for a very massive mode of stellar formation at z ~ 16 (mj n f = 140 M ). These estimates are 
confirmed by results described below. One should note that the outflow is a self-regulating process 
and though it may be large compared with the SFR at a given moment, it will turn off as the 
SFR decreases. Thus we see that the outflow can be ultra-efficient during the early population III 
phase and can still be quite efficient at high redshift even in the normal mode of star formation. 
In contrast, outflows are very inefficient at low redshift. 



2.4. Model parameters 

Before we describe the models in more detail, it will be useful to outline the adjustable param- 
eters that we will consider. While this list of parameters is not exclusive, they are the parameters 
which most affect the results and shape our conclusions. 

1. The minimum mass M m j n of dark matter halos of star-forming structures. As described above, 
the distribution of star-forming structures spans a mass range with some minimal halo mass. 
Structures with masses smaller than M m j n are considered as part of the IGM and do not 
participate in the star formation process. As we will see, when the initial baryon fraction 
is fixed, there is a tight relationship between M m ; n and the redshift at which star formation 
begins. There are few direct constraints on M m j n , and we will consider values from 10 6 - 10 11 
M . 
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2. The initial fraction /b, s truct(^init) of baryons within the structures when star formation begins. 
We will for the most part assume that star formation begins when the baryon fraction in a 
structure is 1%. For a given value of M m j n , eq. (5) can be used to determine the redshift at 
which star formation begins. Thus, we can alternatively specify, Zimtj an d compute the initial 
baryon fraction. 

3. The efficiency of the outflow e. This parameter too, is largely undetermined. Its value can 
be fixed within the model as it sensitively controls the final baryon fraction in structures, 
/b.stmctC- 2 = 0). For models with an exponentially decreasing star formation rate (as we 
consider below), increasing e leads to a decrease in the final value of /instruct (<? = 0) which 
we constrain to lie in the range of 61 ± 18 % (see e.g. Fukugita & Peebles 2004) resulting in 
values of e < 0.02. 

4. The star formation rate ^(t). In principle, there are many possible analytical forms for the 
star formation. For example, it is common to assume that the SFR is proportional to the gas 
mass fraction, a, (or some power of a). As we describe below, we cannot obtain satisfactory 
fits to both the cosmic SFR and the observed type II supernova rates when ^ oc a for the 
normal mode ((0.1 - 100) M Q ) of star formation. In contrast, we find very good fits with an 
exponentially decreasing rate of the form \I> = z/ie - */ r . Similarly, we assume an exponential 
form for the massive population III mode as well, although here we assume a metallicity- 
dependent cutoff so that * = v-ie~ z l ZaVi . The timescale r is determined largely by the 
observed cosmic star formation rate (at low redshift), while the intensity v\ is constrained 
by the observed metal abundances in the ISM and supernova rates, is taken to be as 
large as possible in order to achieve a high ionizing flux at high redshift, while avoiding the 
overproduction of metals which would appear as a prompt initial enrichment of the ISM and 
the IGM. Finally, we assume a critical metallicity of 10~ 4 (Bromm h Loeb 2003; Yoshida 
et al. 2004), though our results do not change significantly if Z cr it were a factor of 10 higher. 

5. The initial mass function of stars <I>(m). As for the SFR, we must specify an IMF for the 
normal and massive modes. In each case, we take <3?(m) oc m~^ 1+x "> for m; n f < m < m sup , 



The normal mode is assumed to contain stars between 0.1 and 100 M . The slope of the IMF 
is constrained by both the type II supernova rate and the metallicity of the ISM. 

6. Specific parameters for type la supernovae. Modelling type la supernovae requires several 
additional parameters which include the mass range [msN ia.min^SN ia,max] of the progen- 
itors, the fraction /sn ia of these progenitors which will produce a type la supernova, and 
time delay AigN i a between the formation of the white dwarf and the explosion (so that for a 
progenitor of mass m, the total delay between the star formation and the type Ia supernova 
is r(m) + At S N ia)- 



with 




(16) 
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3. The normal mode of star formation 

In all of the models considered below, there is a normal mode of star formation comprising 
of stars with masses between 0.1 and 100 M . Once the normal mode is adjusted, a massive 
population III mode is added at high redshift for Z < Z c . Because the massive mode is cut off 
rather early (its main effect is to supply a prompt initial enrichment of the ISM and IGM as well 
as provide an ionizing source at high redshift), most observable constraints are tied to the normal 
mode. Thus, we begin with a more detailed description of the normal mode and its observational 
consequences. 

We listed the main parameters of the model in the previous section. The model consists of a 
superposition of a normal mode of star formation and an early massive mode at high redshift, when 
the global metallicity in the star- forming structures is still very low (population III stars). We will 
describe this mode in section 5. We focus here on the normal mode of star formation. We fix the 
mass range of star formation to mi n f = 0.1 M© and m sup = 100 M Q so that the only parameter 
needed to define the IMF is its slope, namely x±, which is usually estimated to be in the range 1.30 
(Salpeter 1955) to 1.7 (Scalo 1986). 

As noted above, we employ an exponentially decreasing SFR (which is representative of el- 
liptical galaxies) as parametrizations such as a Schmidt-law yield significantly poorer fits to the 
observational data. Hence we take, 

tf(t) = ^iexp(-(t-t init )/ri) , 

where t- m i t is the initial age of the Universe at z- m i t where the star formation starts in the model, t\ 
is a timescale of the order of 2 — 3 Gyr and v\ = /instruct (t)/ T i with f\ being a fraction governing 
the efficiency of the star formation. 

Since many of the observational constraints used to fix the parameters pertain to relatively 
low redshift, we performed a detailed scan of the parameter space including only the normal mode. 
The parameter grid was chosen to be: M min = 10 6 , 10 7 , 10 s , 10 9 , and 10 11 M Q ; e = 0, 10~ 3 , 
2 x 10~ 3 , 3 x 10~ 3 , 5 x 10" 3 , 7 x 10~ 3 , lO" 2 , 2 x 10~ 2 and 3 x 10~ 2 ; v x = 0.1, 0.2, 0.3, 0.4, 0.5, 
0.6 and 0.7 Gyr" 1 ; n = 2.2, 2.4, 2.6, 2.8, 3.0, 3.2 Gyr; x x = 1.3, 1.35 and 1.4. In all cases, we have 
assumed that the onset of star formation begins when the baryon fraction in structures is 1%. We 
have tested that this initial fraction has a very weak impact on the results concerning the normal 
mode of star formation. It is on the other hand of great importance for the population III stars and 
its effect will be studied in section 5. One should also note that there are additional parameters 
associated with the rate of type la supernovae. These are directly constrained by observations and 
will be fixed in section 4. 

To determine the parameters which best fit the observations, we performed \ 2 analysis over 
the parameter grid. Included in the \ 2 analysis are six sets of observational data: (1) The observed 
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cosmic star formation rate up to z ~ 5 (Hopkins 2004). The data were binned and averaged in 
redshift leading to somewhat larger observational uncertainties at a given redshift than typically 
reported for a single measurement; (2) The observed rate of type II supernovae up to z ~ 0.7 
(Dahlen et al. 2004); (3) The present fraction of baryons in structures, /b, s truct (z = 0) ~ 61 ± 18 % 
(Fukugita Sz Peebles 2004); (4) The present fraction of baryons in stars, fb,*(z = 0) ~ 6 ± 6 % 
(Fukugita & Peebles 2004); (5,6) The evolution of the metal content in the ISM and IGM from 
z ~ 5 to z = 0. Specifically, for the ISM we took log 10 Z/Z & = -0.3 ± 0.2 at z = and -0.9 ± 0.2 
at z = 4. For the IGM we took, log 10 Z/Z & = -2.5 ± 0.2 at z = and -3.1 ± 0.2 at z = 4 
(Prochaska et al. 2003; Ledoux et al. 2003; Songaila 2001; Schaye et al. 2003; Aguirre et al. 2004). 
These constraints are chosen so that the metallicity in the ISM (IGM) is always larger (smaller) 
than in the most (least) metallic DLA. The dominant effect of including the massive mode will fall 
on the metallicity of the IGM and will be discussed in section 5 below. 

The results of the x 2 minimization are listed in Table 1 separately for each value of M m j n . 
Having fixed the onset of star formation to correspond to an initial baryon fraction of 1%, each 
value of M m i n leads to a distinct redshift for initial star formation. This value of Zi n n is also given 
in the table. In each case, we found a best fit IMF slope of x\ = 1.3. It should be remembered 
that when we add in the massive mode, these 'best' fits have to be modified in order to take into 
account the metals produced by population III stars. This will affect only the efficiency of the 
outflow, e, as the other parameters control the later phases of evolution. Of the sampled values of 
-^min, ^min = 10 7 M was found to give the best fit. However if the observed SFR at high redshift 
(z > 3) has been underestimated, the best fit model may be at lower values of M m j n . For example, 
we find that M m j n = 10 6 M becomes a better fit if the high redshift SFR data are increased by a 
factor ~ 3. Also shown in the table are the output values for the baryon fraction in structures and 
in stars. Note that none of the models accurately reproduce the stellar baryon fraction at z = 0. 
In each case, we over-produce stars at a level of about la. 

The results of our best fit models are plotted in Figure 3 for each value M m \ n as indicated in 
the upper panel. Because we have fixed the initial baryon fraction in structures, each value of M m - m 
corresponds to a different initial redshift. As one can see, each of the models gives a satisfactory 
fit to the global SFR, save perhaps the case with M m \ D = 10 11 M Q . The data shown (taken from 
Hopkins (2004) has already been corrected for extinction. If we compare the cosmic SFR calculated 
by Croton et al. (2005a) (see also Springel & Hernquist (2003)), we find good agreement at high 
redshift, particularly with the M m \ n = 10 7 M model. However, at z < 5 their SFR is somewhat 
lower than ours by a factor of about 2. Moreover, in our model, about 60 per cent of all stars 
formed by z = 3 compared with 50 per cent at the same redshift in their model. The predicted 
type II supernova rates are also an excellent fit to the existing data at low redshift. The predicted 
rates are substantially higher at higher redshift and lower M m i n . 

As noted above, while the baryon fraction in structures is fit quite well, the baryon fraction in 
stars at z = (or equivalently O?) is somewhat high. This output is clearly linked to the SFR which 
is in fact well fit and could be even higher than shown here. For example, we show in Figure 4 
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Fig. 3. — The cosmic star formation rate, type II supernova rate, baryon evolution and global 
metallicity for the normal mode of star formation (Model 0). Plotted are the results of our best fit 
models for each choice of M m i n . In the top panel, we show the SFR for each model. The observed 
SFR up to z ~ 5 is taken from Hopkins (2004). In the second panel, we show the predicted rate 
for type II supernovae. The observed rates up to z ~ 0.7 is taken from Dahlen ct al. (2004). In 
the third panel we show the evolution of the baryon fraction in structures and in stars (inset). The 
local values at z = are taken from Fukugita <fe Peebles (2004). In the bottom panel, we show 
the evolution of the metallicity in both the ISM and IGM (dashed curves). The metallicity of 100 
Damped Lyman-a systems as measured by Prochaska et al. (2003) is also plotted. 
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Table 1: The best fit models for the normal mode of star formation (Model 0). Column 1 indicates 
the input value of the minimum mass for star forming structures. Column 2 is derived from column 
1, having assumed that /& = 1% when star formation begins. In columns 3, 4, and 5, we show the 
results of the minimization for each value of M m j n . The slope of the IMF is x = 1.3 for all these 
best models. Among all these models, the best fit is obtained for M m \ D = 10 7 M (bold face). For 
this reason, this model is used as a reference in the figures. Columns 6 and 7 contain the output 
values of the baryon fraction in structures and in stars respectively. 

resulting predictions when the astration rate V\ has been doubled for our best fit model with 
-^min = 10 7 M Q . As one can see, the baryon fraction in stars now rises to nearly 30% and clearly 
disagrees with observational determinations of this quantity. Finally, coming back to Figure 3, we 
also show the evolution of the metallicity in both the ISM and IGM (dashed curves). While the ISM 
metallicity rises very quickly initially, evolution is more gradual in the IGM where it is controlled 
by the outflow. Furthermore, we see that the evolution of the metallicity in the IGM is dependent 
on M m i n . Since models with large M m j n begin rather late (e.g. the model with M m j n = 10 11 M Q 
begins at z ~ 6.5) these models have difficulty in producing a suitably large metallicity in the IGM 
at redshift 2-3. Note that we have only included the normal mode here, and contributions from 
population III stars will only serve to increase the metallicity in both the ISM and IGM. 

4. Supernova rates 

Among the important constraints imposed on these models of cosmic chemical evolution are 
the rates of type la and II supernovae. These are intimately linked to the choice of an IMF and 
SFR and represent an independent probe of the star forming universe at high redshift. While type 
II rates directly trace the star formation history, there is a model dependent time delay between 
the formation (and lifetime) of the progenitor star and the type la explosion. Since existing data 
is available only for relatively low redshifts, we discuss the predicted SN rates in the context of 
Model (the normal mode of star formation). Indeed, at these redshifts the massive mode is not 
operational. 
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Fig. 4. — As in Figure 3, we show the best fit model with M m i n = 10 7 M (thin line) and the same 
model where the astration rate v\ has been arbitrarily multiplied by a factor 2 (thick line). Note 
the new high value of the corresponding local fraction of baryons in stars and the increase in the 
peak value of the SN rate at a redshift of about 3. 
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4.1. Type II supernovae 

The predicted rates for type II supernovae were discussed briefly in the previous section, where 
we showed the results as a function of redshift for several choices of M m i n in Figure 3. In Figure 5, 
we focus on the low redshift range and concentrate on the best fit model with M m \ n = 10 7 M Q . 
The data shown are taken from the Great Observatories Origins Deep Survey (GOODS) (Dahlen 
et al. 2004, and references therein). Of the 42 supernovae observed, 25 are associated with type la 
and 17 with core collapse supernovae (type II as well as type Ib,c). 

The GOODS data for core collapse supernovae have been placed in two bins at z = 0.3 ± 0.2 
and z = 0.7 ± 0.2. Their results (which have been corrected for the effects of extinction) show SN 
rates which are significantly higher than the local rate (at z = 0) which is taken from Cappellaro 
et al. (1999). Since the timescale between star formation and the core collapse explosion is very 
short (particularly for very massive stars) , there will be no contribution of Population III stars to 
type II rates at z < 1. 

As remarked above, the SN type II rate is directly related to the overall SFR and therefore 
to the assumed astration rate v\ and our best fit model does an excellent job at describing the 
existing data. Note that, had we chosen a Schmidt law (\l/ cx a gas ), it would not be possible to 
simultaneously fit both the local [z = 0) data as well as the GOODS data. In general, models 
with a Schmidt law can not account for a large enhancement in the SN rate. Therefore once the 
peak of the SFR has been fixed (at z « 1), the local SN rate is too large. The slope of the type II 
rate (with respect to redshift) is also dependent on the IMF slope, x\. In Figure 6, we show this 
dependence displaying the predicted type II rates for x\ = 1.1, 1.3 and 1.5. Of course the choice of 
v\ and x\ also has a strong effect on the metallicity and detailed chemical history as discussed in 
more detail in the next section when Population III stars are included in the analysis. 

It is important to note that while the data at z < 1 is well described by the model, the same 
model predicts much higher SN rates at higher redshift. Indeed, the SN type II rate peaks at a 
redshift z « 3 at a rate which is nearly 8 times the observed rate at z « 0.7. It will be interesting to 
see whether future data will be able to probe the SN rate at higher redshift. These elevated rates 
may be detectable in spite of the expected increase in dust extinction due to the early production of 
metals. An alternative to the direct detection of SN at high redshift is the indirect detection of SN 
through neutrinos. The predicted neutrino background is very sensitive to the assumed chemical 
evolution model (Ando k, Sato 2004; Iocco et al. 2005; Daigne et al. 2005). 

4.2. Type la supernovae 

As in the case for type II supernovae, the data (also taken from GOODS (Dahlen et al. 2004)) 
is binned into four redshift bins at z = 0.4, 0.8, 1.2, and 1.6, each with a spread of ±0.2. These are 
shown in the lower panel of Figure 5 along with data from Reiss (2000), Hardin et al. (2000), Pain 
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Fig. 5. — The evolution of the ratio of type la to type II supernova event rates is shown for the 
best fit model with M m i n = 10 7 M in the upper panel. The observed local values are taken from 
Cappellaro et al. (1999). The type II SN rate as a function of redshift is shown in the middle panel. 
The observed rates at z ~ 0.3 and z ~ 0.7 are taken from Dahlen et al. (2004). In the bottom 
panel, we show the type la SN rate as a function of redshift. The observed rate up to z ~ 1.6 is 
taken from Reiss (2000), Hardin et al. (2000), Pain et al. (2002), Madgwick et al. (2003), Tonry 
et al. (2003), Strolger et al. (2004), Blanc et al. (2004) and Dahlen et al. (2004). 
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Fig. 6. — The effect of the IMF slope on the SNII rate as a function of redshift. The reference 
model (plotted as the thin curve) assumes M m \ n = 10 M and an IMF slope of x\ = 1.3). For 
comparison, two other models with slopes x = 1.1 and x = 1.5 are also plotted. The data shown 
are the same as in Figure 5. 
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et al. (2002), Madgwick et al. (2003), Tonry et al. (2003), Strolger et al. (2004), and Blanc et al. 
(2004). Once again, the local (z = 0) data is taken from Cappellaro et al. (1999). Unlike the type 
II data, there appears to be a definite peak in the rate of type la supernovae around z ~ 1. 

As the progenitors of type la supernovae are stars with intermediate masses between 2 and 
8 M , the rate for SN type la is controlled entirely by parameters within the normal mode of 
star formation. There are however two additional parameters that must be adjusted to fit the 
observed data, namely the fraction of intermediate mass stars that explode as SN type la, and the 
time delay between the death of the progenitor and the SN explosion. In most previous models 
of Galactic chemical evolution, this delay was taken to be of order 1 Gyr. Recently, many studies 
have been devoted to the progenitors of SNIa supernovae, including the analysis of time delays and 
the mass ranges (Han & Podsiadlowski (2004), Yungelson (2004), Garnavich & Gallagher (2005) 
and Belczynski et al. (2005)). The GOODS data are now consistent with significantly larger delays 
(of order 4 Gyr, including the lifetime of the star) and exclude delays less than 2 Gyr at the 95 % 
confidence level (Strolger et al. 2004; Dahlen et al. 2004). In the context of the normal mode (Model 
0), we find that these observations are well reproduced if ~ 1 % of intermediate mass stars lead 
to type la supernovae and if the typical delay between the formation of the white dwarf and the 
explosion is ~ 3 — 3.5 Gyr. The resulting comparison between our theoretical prediction and the 
data is shown in the lower panel of Figure 5 where the effect of the time delay is quite apparent and 
is chiefly responsible for the peak in the SN rate at z ~ 1. While this comparison is made using 
our best fit model with M m j n = 10 7 M , the resulting rate for type la supernovae would be quite 
similar for other choices of M m j n . For an alternative approach see Scannapieco & Bildsten (2005); 
Mannucci et al. (2005). 

To test the sensitivity of the result to the time delay, we show in Figure 7 the predicted SNIa 
rates using a time delay of 2.0 and 4.0 Gyr, compared with our best fit value of 3.2 Gyr. Also 
shown in Figure 7 (upper panel) is the sensitivity to the adopted intermediate mass range. Our 
best fit model assumes a range of 2 - 8 M . Also plotted are models where the lower end of the 
mass range is 1.5 and 3 M . 

Finally, in the upper panel of Figure 5, we compare our predictions of type la to type II super- 
novae. We see that these are quite consistent with the local rate (Cappellaro et al. 1999). The ratio 
drops off at redshifts z > 1 due to the time delay in producing type la events. Furthermore, when 
SN la explosions do occur, structures are larger and mass outflows are suppressed (see discussion 
in the next section). Thus, as a consequence of the SN la time delay, we predict that while ~ 50 % 
of iron in structures is produced by type la supernovae, this fraction is only ~ 10 % in the IGM. 



- 21 - 



in 

2 25 



20 
15 



13.7 



25 - 



20 



15 



10 - 



5 - 







10 - 



5 - 



- 



Time (Gyr) 
5 

I 



2.5 



m 



SN Ia.min 



= 1.5 2.0 3.0 




At SNIa = 40 3 2 20 




t± I I I I I I I I I I I I I I I I I I I I I I I I L 







0.5 



1 1.5 
Redshift 



2.5 



Fig. 7. — The effect of the progenitor mass range and of the time delay on the SNIa rate. The thin 
curves correspond to the reference model with M m - m = 10 7 M , a minimal progenitor mass of 2 
Mq, and a time delay of 3.2 Gyr. The sensitivity to the lower end of the intermediate mass range is 
shown in the upper panel where two other models with a minimum mass of 1.5 and 3 M are also 
indicated. The sensitivity to the time delay is shown in the lower panel where two other models 
with time delays of 2 and 4 Gyr are also indicated. The data shown are the same as in Figure 5. 



- 22 - 



5. Population III stars 



5.1. Necessity of an initial starburst of massive stars 



As shown in Daigne et al. (2004), the normal mode of star formation, labelled Model there as 
well as here, is not capable of reionizing the Universe at high redshift. Indeed, at a redshift z = 17, 
only 1.6 ionizing photons per baryon are available compared to the requisite value of approximately 
20 assuming a dumpiness factor, Ch n = 10 (Ricotti & Ostriker 2004). It was also shown in Daigne 
et al. (2004) that Model alone is not capable of explaining the observed abundance patterns in the 
extremely iron-poor stars, CS 22949-037 (Depagne et al. 2002; Israelian et al. 2004), HE 0107-5240 
(Christlieb et al. 2004; Bessell et al. 2004), HE 1327-2326 (Frebel et al. (2005)) and G 77-61 (Plez 
& Cohen (2005)). 

The necessity of a massive mode of star formation, active at high redshift has become an 
integral part of our emerging picture of the growth of galactic structures. However, there is an 
active debate as to the exact nature of the massive mode. As in Daigne et al. (2004), we consider 
three possibilities for the massive mode. 1) stars with masses in the range 40 - 100 M©, Model 
1; these stars terminate as type II supernovae. 2) stars with masses in the range 140 - 260 M , 
Model 2a; these stars terminate as pair- instability supernovae (PISN). 3) stars with masses in the 
range 270 - 500 M , Model 2b; these stars terminate as black holes through total collapse and do 
not contribute any metal enrichment. 

In all cases, we assume a bimodal birthrate function of the form 



where the two IMFs are normalized independently as in eq. (16). The SFR, vl^ is now expressed 
as a function of metallicity and cuts off once a critical metallicity is reached, as described below. 
Furthermore, we assume the same slope of the IMF, x\ = X2 = 1-3 so that $i and $2 differ 
only in their mass range. In principle, we can choose to start the normal mode of star formation 
either simultaneously with the massive mode, or sequentially, that is, when the Z > Z C rit- In the 
latter case, one can argue that the initial injection of metals by Pop III stars is responsible for 
the formation of the first extremely metal poor Pop II (or Pop II. 5) stars (Mackey et al. 2003; 
Salvaterra et al. 2004; Johnson & Bromm 2005). Since critical metallicity is achieved very rapidly 
(within 3 Myr), our results for these two choices are almost indistinguishable. 

As the cooling process of the gas depends strongly on its chemical composition, it is believed 
that the evolution of the mass range of the IMF is mainly governed by the global metallicity (Fang 
& Cen 2004). As noted above, we assume a transition from population III to the normal formation 
mode at a critical metallicity Z a \ t (Bromm h Loeb 2003; Yoshida et al. 2004) by defining the SFR 
of the massive mode by 



B(m,t,Z) = *i(m)tfi(t) + $ 2 (m)tf 2 (2) 



(17) 



V(t) 2 = v 2 exp(-Z/Z CTit ) , 
with v 2 = instruct (*)• We adopt Z crit /Z Q = 10~ 4 . 



(18) 
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Because the massive stars adopted in Models 1 and 2a are efficient in producing heavy elements, 
the duration of the massive phase in these models is relatively brief. In contrast, no heavy elements 
are produced in Model 2b, and therefore, the massive mode continues to affect the evolution of 
structures until the metallicity of the normal mode reaches the critical value, which is also relatively 
fast, due to the short lifetimes of massive stars. 

Most of the constraints discussed above must be applied at relatively low redshift [z < 5). 
As a result they fix the parameters of Model 0. For the most part these are unchanged, with one 
important exception. The metallicity of the IGM now receives two distinct contributions. The first 
from the massive mode, which will appear as a prompt initial enrichment, and the second from the 
outflows of the normal mode which contribute at lower redshifts. In principle, these contributions 
can be distinguished with precise abundance data as a function of redshift. If the metallicity in 
the IGM originates from the massive mode, we expect IGM abundances which are constant with 
respect to z. On the other hand, if there is a significant contribution from the outflows due to the 
normal mode, we should see abundances which vary with redshift. Unfortunately, the current data 
does not allow us to distinguish between these two possibilities. Indeed, one would expect that 
the IGM contains contributions from both modes, but the relative contribution in fact becomes 
another parameter of our model. 

In what follows, we will consider two possibilities for the role of the massive mode in the 
enrichment of the IGM. In one case, we will assume that nearly all (~ 90%) of the IGM metal- 
licity originates from massive stars. Increasing the contribution of the massive mode allows us to 
maximize its ionizing efficiency In section 3, we discussed the role of the normal mode. There, 
we had effectively set i>2 = 0, and the efficiency of outflow, e, controlled both the baryon fraction 
in structures and the metallicity of the IGM. When v?, is maximized, we will clearly be forced to 
reduce the value of e to avoid the overproduction of metals in the IGM. Despite the low value of 
e, outflow is more efficient at early times and the enrichment of the IGM is also more efficient due 
to the high level of SNe (large v<i). Consequently, in this case, most of the IGM metals comes from 
pop III stars and we expect that Zjqm will be constant with z. Note that although the value of 
Z CT i t is very small, ISM metallicities and even IGM metallicities much greater than Z cr jt can be 
produced by the massive mode. This is because of the finite lifetimes of the massive stars. Although 
these lifetimes are short, they are nevertheless of order a few Myr and hence metals continue to be 
injected into the ISM (and IGM) after Z > Z crit is reached and the formation of new massive stars 
is quenched. 

We will also consider an intermediate case where the IGM receives equal contributions from 
the massive and normal modes at z = 2.5. In this case, v<i is again chosen to maximize the efficiency 
of reionization, and e will be closer to the value found in section 3. 

In Table 2, we display the four parameters for Models 1, 2a, and 2b, for both sets of assumptions 
concerning the contribution of Pop III to the IGM. Models 1, 2a, and 2b refer to the case where 
the IGM receives equal contributions from the normal and massive modes (at z w 2.5), whereas 
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Table 2: Parameter values for the massive starburst Models 1, 2a and 2b. Column 1 indicates the 
model number and column 2 the input value of the minimum mass for star forming structures. In 
columns 3 and 4, we show the outflow efficiency and massive mode SFR. 

Models, le and 2ae, refer to the extreme case where 90% of the metallicity at z ~ 2.5 is derived 
from the massive mode. Note that there is no Model 2be, as no metals are ejected into either the 
ISM or IGM from massive mode stars in this case, and all of the low redshift metallicity is due to 
the normal mode. As noted above, the outflow efficiencies for Model 1 are within a factor of 2 of 
those for Model (cf. Table 1). In contrast, the efficiencies for Model le are significantly lower. 
However as one can see in Table 2, the low outflow efficiency is compensated for by a large SFR 
(1/2) greatly increasing the models capacity for reionization. The final baryon fraction in structures 
and stars is very similar in these models to that found in Model 0. In general, instruct is larger 
than the values in Table 1 by 1 - 2 % for Models 1, 2a, and 2b, and by 2 - 4 % for Models le 
and 2ae. /b * is unchanged from Model 0. As we will see, there is also a close correlation between 
the intensity of the SFR associated with the massive mode, reionization and the metal enrichment 
provided by these Pop III stars. 

Once the ratio of massive/normal mode contribution to the IGM is adopted, we apply the 
observational constraints to fix the values of e and v<i- In addition to the constraint based on 
the overall metallicity of the IGM, we must require that the massive mode produces a sufficient 
number of ionizing photons. We will also look at the evolutionary history of several individual 
element abundances including C, O, Si, and Fe. Finally, we compare the results to the observed 
abundances of several extremely iron-poor stars CS 22949-037 (Depagne et al. 2002; Israelian et al. 
2004), HE 0107-5240 (Christlieb et al. 2004; Bessell et al. 2004), HE1327-2326 (Frebel et al. 2005), 
and G 77-61 (Plez & Cohen 2005). They represent the mean abundance in the Universe at this 
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cpoch. 

5.2. Model 1 

We begin by discussing the results of the bimodal model which combines the normal mode 
(Model 0) with the massive mode denoted as Model 1 which includes stars with masses in the range 
40 -100 M 0. In this model (as well as in Model 2a described below), star formation begins at a 
very high rate and falls precipitously as metals are injected into the ISM. In inhomogeneous models 
of structure formation, we would expect this narrow burst to broaden (Scannapieco et al. 2003). 
As in the case of Model 0, the onset of star formation is determined by M m - m , and the initial value 
for the baryon fraction in structures (fixed to be 1 %). Results for the SFR, baryon fraction and 
metallicity are shown in Figure 8 for the intermediate case where Model 1 stars contribute 50% to 
the IGM metallicity at z = 2 to 3 and 20% at z = 0. With exception of the metallicity, shown in the 
third panel, the effect of the Population III stars is minor as the SFR and baryon fraction is very 
similar here to the case of Model shown in Figure 3. Note that at high redshift, the fraction of the 
baryons in massive stars is about 10 -3 as seen in the insert to the middle panel of Figure 8. As one 
can see in the lower panel, the metallicity in the ISM reaches values far in excess of Z cr [ t due to the 
finite lifetime of the massive stars relative to the speed at which the metallicity is attained. Notice 
also that once the metallicity from Pop III stars is produced, the ISM metallicity later decreases 
as a result of the accretion of metal-free gas as the structures grow. 

The rate of accretion as well as the rate of outflow is shown in Figure 9. We see here that 
outflows from Pop III stars are very efficient at high redshift but their duration is very short due 
to the rapid increase in metallicity in the ISM. It is important to remember that the mass flux 
of outflows is very sensitive to e. As a consequence of the luminosity density obtained by the 
Sloan Digital Sky Survey (SDSS, Blanton et al. (2003)) and by the 2dF galaxy redshift Survey 
(Croton et al. 2005b) together with the evaluation of the mass rate triggered by galactic winds 
(Veilleux et al. 2005) we estimate that the actual mass flux of outflows is of the order of 0.01 to 0.1 
M0/yr/Mpc 3 . At z = 2, the results for Model 1 show outflows at the lower end of this range. 

Among the chief motivating factors in developing a model of cosmic chemical evolution is the 
early reionization of the Universe. Several studies suggest that an early burst of star formation 
as in Model 1, is sufficient (Venkatesan et al. 2003; Wyithe & Loeb 2003b; Venkatesan & Truran 
2003; Tumlinson et al. 2004; Daigne et al. 2004). In Figure 10, we show the number of ionizing 
photons per baryon produced for each of our choices of M m [ n . The procedure for computing this 
stellar ionizing flux is explained in DOVSA. It is important to remember that only a fraction 
/esc of these UV photons will escape the structures and therefore be available to ionize the IGM. 
The effective value of / e sc is poorly known but could vary from about 1 to 30 %. The minimum 
number of photons required for complete reionization (see DOVSA) is also plotted for three possible 
dumpiness factors. The ionizing potential clearly increases with M m - m and decreasing redshift. The 
ratio of this minimum number (dashed line) to that for the stellar ionizing photons (solid line) 
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Fig. 8. — As in Figure 3 but with an added initial starburst with stellar mass in the range 40-100 
Mq (Model 1). See Table 2 for the parameters used. 
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Fig. 9. — The rate of structure formation and outflow in Model 1. The mass flux corresponding to 
the structure formation is given by the solid curves and the rate of outflow from the structures to 
the IGM is shown by the dotted curves as a function of redshift for all models shown in Figure 8. 
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gives the minimum fraction / esc necessary to fully reionize the IGM. It appears that all models 
are able to fully reionize the IGM for a dumpiness factor Ch ii ~ 10 and / esc ~ 25% — > 1% for 
M m i n = 10 6 — > 10 11 Mq. When the dumpiness factor increases, only models with high M m i n 
are still able to fully reionize the IGM with / csc < 30%. Typically, if Ch ii ~ 30 (100), we need 
M m i n > 10 7 M (10 9 M Q ). There is therefore a tendency to favor a high M min and low reionization 
redshift if the dumpiness factor is high. 

In Figure 11, we show the dependence of these results on the choice of Z c , the critical metal- 
licity which shuts down massive star formation, and the initial baryon fraction in structures which 
allows star formation. In all of the models discussed up to now, we have taken Z c = 1O _4 Z and 
/b,struct(^init) = 1%- We see in Figure 11 that results are almost independent of Z c . Increasing 
Z c leads to a modest change in the number of ionizing photons per baryon. The reason is that 
the first generation of stars formed are in fact capable of producing a metallicity much larger 
than either 10~ 4 or 5 x 10~ 3 Z Q , and massive star formation is rapidly terminated once these stars 
pollute the ISM. On the other hand, the results for the number of ionizing photons per baryon 
are quite sensitive to the choice of instruct (zinit) as is the initial redshift for star formation. Our 
standard choice for /b struct ( 2 imt) of 1% corresponds to our conservative estimate of the minimum 
baryon fraction where sufficient dissipation occurs to allow star formation. Most galaxy formation 
simulations adopt a value of 10% (cf. (Abadi et al. 2003)). 

Decreasing /instruct (zinit) by a factor of 10 to 0.1% allows star formation to begin earlier (in this 
case at a redshift z ~ 21 for M m \ D = 10 7 M ), when the IGM is still quite dense, and reionization 
does not occur, even if Ch ii = 10. If /b,struct(zinit) is increased to 5%, star formation occurs later 
(z ~ 11 for M m i n = 10 7 M Q ), when the IGM is less dense, and reionization occurs quite easily, even 
if Ch ii = 100 (only a fraction / csc ~ 13% is needed in this case). 

Next, we show results for the evolution of C, O, Si, and Fe as a function of redshift in Model 
1. These are found in Figures 12 and 13 (the yields vary with metallicity and come from Woosley 
& Weaver (1995)). Overall, one sees that the chemical abundances are well reproduced. In the 
cases of O and C, we see quite clearly the effects of the hierarchical growth of structure. After the 
initial burst of star formation by massive stars, the ISM abundances within structures is diluted 
as metal-free baryons are accreted. As the normal mode of star formation begins to eject metals, 
these abundances subsequently begin to rise again. Unfortunately, observations at sufficiently high 
redshift do not exist at present to test this predicted feature. Because of the rapid ejection of 
metals in Model 1, we are able to achieve the high abundances of C and O seen in the extremely 
iron-poor stars indicated in Figure 12, in contrast to what is possible in Model 0, especially for C 
(see DOVSA). 

Model 1 predictions in the IGM agree with the estimate of the oxygen abundance in the Lyman 
alpha forest (Telfer et al. (2002), Simcoe et al. (2004), Aracil et al. (2004), Bergeron & Herbert-Fort 
(2005)). Note that only OVI is directly observed. Therefore, the estimate of the total O abundance 
is limited by the uncertainties of the ionization state of this element. We can reasonably consider 
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Fig. 10. — The cumulative number of ionizing photons emitted by stars is plotted as a function of 
redshift for all models of Figure 8. The minimum number of photons per baryon necessary to fully 
reionize the IGM is plotted as a function of redshift by thin dashed lines for three different values 
of the dumpiness factor Ch ii = 10, 30 and 100. 
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Fig. 11. — As in Figure 10 for M m - m = 10 7 M©, we plot the number of ionizing photons per 
baryon. The dotted curve corresponds to an increase of Z c to 5 x IO^^Zq, ie. by a factor of 50. 
The dashed curves show the number of photons per baryon for two different choices for the initial 
baryon fraction in structures 0.1 % or 5 % as compared with the nominal choice of 1 %. 
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Fig. 12. — The evolution of the Oxygen and Carbon abundances in Model 1. The mass fraction of 
Oxygen (top) and Carbon (bottom) as a function of redshift is plotted both in the ISM of structures 
(solid curves) and in the IGM (dashed curves) for all models of Figure 8. The observed abundances 
are taken from Pettini et al. (2002) (filled triangles: Oxygen abundance in 10 DLAs), Shapley 
et al. (2004) (open boxes: Oxygen abundance in 7 massive galaxies), Simcoe et al. (2004) (Big 
dot: Oxygen VI in IGM. The arrow indicates the estimated value of the total Oxygen abundance), 
Songaila (2001) (circle: Carbon IV abundance in the IGM) and Aguirre et al. (2004) (big dot: total 
Carbon abundance in the IGM). Horizontal thin dashed lines indicate the measured abundances in 
the following four very metal-poor halo stars: CS 22949-037 (Oxygen and Carbon, from Depagne 
et al. (2002)), HE 0107-2240 (Oxygen and Carbon, from Bessell et al. (2004)), HE 1327-2326 
(Oxygen and Carbon, from Frebel et al. (2005)) and G 77-61 (Carbon, from Plez & Cohen (2005)). 
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Fig. 13. — As in Figure 12, the evolution of the Iron and Silicon abundances in Model 1. Observed 
abundances are taken from Prochaska et al. (2003) (crosses: Iron and Silicon abundances in 100 
DLAs), Ledoux et al. (2003) (triangles: Iron and Silicon abundances in 33 DLAs and sub-DLAs), 
Songaila (2001) (circle: Silicon IV abundance in the IGM) and Aguirre et al. (2004) (big dot: total 
Silicon abundance in the IGM). Horizontal thin dashed lines indicate the measured abundances in 
the following four very metal-poor halo stars: CS 22949-037 (Iron and Silicon, from Depagne et al. 
(2002)), HE 0107-2240 (Iron and Silicon, from Bessell et al. (2004)), HE 1327-2326 (Iron, from 
Frebel et al. (2005)) and G 77-61 (Iron, from Plez & Cohen (2005)). 
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that [0/H] can be as high as -2.5 (P. Petitjean, private communication). OVI was recently used 
to ascertain the distribution of metals in the IGM (Tumlinson & Fang 2005). Similarly, the IGM 
abundances of C agree well with the determination by Aguirre et al. (2004) and are (as expected) 
in excess of the CIV determinations of Songaila (2001, 2005). In contrast the models currently fit 
the SilV abundances and underproduce the total Si abundance based on the adopted yields of type 
II supernovae in Model 1. Indeed, we see from Fig. 13 that the Si abundance derived from SilV 
observations (Aguirre et al. 2004) is clearly underproduced by stars in the mass range defined by 
Model 1. As we will see below in section 5.5, this translates in to a predicted value for [Si/C] which 
is low compared with the determination of Aguirre et al. (2004). Furthermore, we note that as 
M m i n is increased, it becomes more difficult to reproduce the IGM abundances. 

To test the importance of the normal mode of star formation on the evolution of the chemical 
abundances, we also consider a set of models in which the efficiency of outflow is reduced, thereby 
increasing the relative contribution of Population III stars to the overall metallicity in the IGM. To 
compensate for the reduction in outflow, the SFR of the massive mode is increased (cf. Model le in 
Table 2) which enhances the ability of the population to reionize the IGM. In Figure 14, we show the 
SFR, baryon fraction, and metallicity for Model le compared with the analogous result for Model 
1, restricting attention to M m \ n = 10 M . With the exception of e, other Model parameters 
are left unchanged. The difference in the SFR between Models 1 and le is predominantly in the 
initial burst of massive star formation. The massive SFR is typically 5 times larger initially, while 
at lower redshifts the SFR, which is determined by v\ is identical. Similarly, the baryon fraction 
in stars is the same for both Models 1 and le, but le exhibits a strong peak (close to = 1), 
which later relaxes due to hierarchical growth. Unfortunately it is not possible to observationally 
test for these differences. In contrast, however, the overall metallicity does present us with a real 
test (requiring suitably accurate data) between Models 1 and le. In Model 1, the IGM metallicity 
is clearly changing with redshift. For example, when M m - m = 10 7 M Q , the metallicity increases by 
a factor of about 20 from a redshift of 16 to and by a factor of about 2 from a redshift of 5 to 
0. In Model le, the evolution of Z is very nearly constant in the IGM, with its final value fixed at 
very high redshift. 

One of the key benefits of Model le is its ionizing potential relative to Model 1. In Figure 15, 
we compare the number of ionizing photons per baryon in Models 1 and le. As one can see, 
whereas Model 1 (with M m - m = 10 7 M ) provided the minimum number of ionizing photons when 
Ch ii = 30, Model le produces a factor of 5 times more ionizing photons per baryon at high redshift. 

As explained earlier, to achieve a high ionization efficiency, we increased the SFR parameter 
V2- However, in order to avoid overproduction of IGM metals formed by Pop III stars, we are 
required to adjust e downwards. Consequently, while outflow is very efficient at very high redshift, 
it is suppressed at later times due to hierarchical growth and is unable to eject metals coming from 
subsequent structures. The mass flux of outflows in Models 1 and le are compared in Figure 15. 
Recall that observations indicate a flux of 0.01 - 0.1 M /yr/Mpc 3 . We see that the predicted 
outflow in Model 1 is ten times higher than that in Model le. Model 1 seems to be in better 
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Fig. 14. — The star formation rate, baryon evolution and global metallicity for Models 1 (thin 
lines) and le (thick lines and labelled extr.) with M Ui \ a = 10 7 M . 
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Fig. 15. — The reionization potential and outflow rate for Models 1 (thin lines) and le (thick lines 
and labelled extr.) with M min = 10 7 M . 
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agreement with the data than Model le, which is probably insufficient to fit the high outflows 
observed in galaxies at z = 2 (see also Bertone et al. (2005); Rupke et al. (2005)). 

Despite the advantage in ionization potential, it is also difficult to reconcile the chemical history 
of Model le with IGM observations. In Figures 16 and 17, we compare the evolution of C, O, Si, 
and Fe in Models 1 and le. The early burst of star formation in Model le produces a prompt initial 
enrichment of the IGM in C and O, which show very little evolution at lower redshifts. In contrast, 
the evolution of Si and Fe is particularly poor, as this models fails to reproduce the observed IGM 
abundances of these elements at low redshift. 



5.3. Model 2a 

Next we consider a massive mode made of stars with masses in the range 140 - 260 M© which 
explode as pair-instability supernovae. In Figure 18, we compare the star formation rate, ionization 
potential and metallicity in Model 2a with that of Model 1. In Model 2a, the initial rapid burst 
of star formation leads to large metallicities at high redshift. In the ISM, these are diluted by 
hierarchical growth. In the IGM, where there is no dilution, the metallicity remains relatively 
flat. Our results confirm that PISN models are less efficient at reionization than a more standard 
IMF as in Model 1. Indeed, even for a dumpiness factor, Ch ii ~ 10, Model 2a would require a 
photon escape fraction in excess of 50% placing this model in a difficult position to explain the 
early reionization of the IGM. 

We have also considered a more extreme version of Model 2a (2ae) in which approximately 
90% of the IGM oxygen abundance at z = 2.5 originates in Pop III. As with the relation between 
Models 1 and le, Model 2ae has a star formation rate characterized by vi which is about 5 times 
larger than Model 2a. This enhancement is seen in the early SFR is shown in Figure 19 where the 
evolution of the baryon density and overall metallicity are also compared with Model 2a. 

In the more extreme case, (Model 2ae), the number of ionizing photons per baryon (shown in 
Figure 20) is increased to about 100 at z ~ 16, making this version of the PISN model acceptable 
based on reionization. However, as was the case in Model le, the mass flux of outflows at low 
redshift is very small in Model 2ae. 

The evolution of C, O, Si, and Fe are shown in Figures 21 and 22. This model is acceptable 
for all nuclei (see however section 5.5 for a comparison of abundance ratios with the those observed 
in iron-poor stars), and predicts little evolution in the IGM. 

5.4. Model 2b 

Finally, we turn to the case of Model 2b, where the massive mode is made of stars with 
masses in the range 270 - 500 M which collapse directly to black holes without any ejection of 
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Fig. 16. — The evolution of the O and C abundances in Models 1 (thin lines) and le (thick lines 
and labelled extr.) with M min = 10 7 M . 



-38- 



13.7 5 2.5 



Time (Gyr) 
1 0.5 



0.21 



G 
o 

• t— I 

-!-> 

o 

CO 

u 

«4-H 

w 
w 
cC 
S 



10- 3 

10" 4 

io- 5 

Et, 

S 10" 6 

io- 7 

10" 8 

io- 9 

IO" 3 

io- 4 
io- 5 



£ io- 6 

CM 



io- 7 

IO" 8 

io- 9 



100 DLA systems (Prochaska et al. 2003) J 
49 DLA systems (Ledoux et al. 2003) 



extr. 




CS 22949-037 
G 77-61 
HE 1327-2326 



IGM (Songaila 2001) □ 
21 DLA systems (Ledoux et al. 2003) 

Aguirre et al. (2004) • 

extr. 



CS 22949-037 




10 

Redshift 



Fig. 17. — The evolution of the Si and Fe abundances in Model 1 (thin lines) and le (thick lines 
and labelled extr.) with M min = 10 7 M . 
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Fig. 18. — The star formation rate, reionization potential, and global metallicity for Models 1, 2a, 
and 2b with M min = 10 7 M Q . 
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Fig. 20.— 



As in Figure 15 for Model 2a. 
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As in Figure 16 for Model 2a. 
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heavy elements into either the ISM or the IGM. The star formation rate, ionization potential and 
metallicity in Model 2b are also shown in Figure 18. As one can see, the SFR and the number of 
ionizing photons per baryon are quite high relative to Models 1 and 2a. Because these stars do not 
produce any metals, the chemical evolution of this model is identical to that of Model 0. This is 
mostly problematic in regard to the very iron-poor stars whose abundances could not be explained 
in this case. 

5.5. Abundance Ratios 

Before concluding, we examine and compare a number of abundance ratios in Models 1 and 
2a. In particular, we consider the evolution of O/Fe, C/Fe, and Si/C. These ratios are of particular 
interest when trying to model the very iron-poor stars which exhibit anomalously high ratios of 
O/Fe, C/Fe and to some extent Si/C. The latter is of potential interest at lower redshifts where 
some derivations of the IGM abundance of Si/C are relatively high coming from SilV and CIV 
observations. 

As shown in Figure 23, the extremely iron-poor stars observed over the last several years show 
values of [O/Fe] between 2-4, which is significantly higher than typical Population II values which 
range between 0.5 and 1. A similar pattern is seen for [C/Fe] which is also very large in these 
stars. Model 1 reproduces these patterns quite well because high C,0/Fe ratios are expected in 
massive type II supernova explosions. The Si/C is also suppressed, and the fit to HE 0107-5240 is 
also quite good whereas the model somewhat underpredicts the Si/C ratio for CS 22949-037. This 
result reinforces the notion that these stars are indeed very old and reveal the composition of the 
gas at the end of Population III. In contrast, Model 2a and PISN are not capable of achieving the 
element ratios observed in these stars. 

At low redshift, one can attempt to use the IGM abundances of Si/C to try to differentiate 
between models. From Figure 23, we see that Model 1 predicts a [Si/C] ratio of order -0.25 whereas 
in Model 2 it is closer to 1 at a redshift between 2 and 4. This is due to the large Si yields found 
in PISN models (Heger & Woosley 2002). There have been some recent observations of ionized Si 
and C (Songaila 2001; Schaye et al. 2003; Aguirre et al. 2004) from which a model-dependent ratio 
was inferred (Haardt & Madau 1996). Indeed the result is highly dependent on the assumption 
of the UV background model. The most extreme case is that of a UV background powered solely 
quasars. In this case, [Si/C] = 1.48, far exceeding the model predictions found here. When a 
quasar plus galaxy model is considered, the ratio drops to 0.77, which is close to the result of 
Model 2a. When the UV flux is softened at high redshift, the ratio drops to 0.48. This ratio may 
present a serious challenge to Model 1 which clearly underproduces Si. We note however, that the 
resulting Si abundance is sensitive to several parameter choices. For example, lowering m in f to 
20 M results in more Si (but aversely affects the C,0/Fe ratios as shown in DOVSA. Increasing 
the slope of the IMF and more importantly decreasing the massive mode astration factor, zv 2 , both 
lead to enhanced Si production. The latter effect can be seen in Fig. 17, where we show the effect 
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Fig. 23.— The evolution of the [O/Fe], [C/Fe] and [Si/C] ratios in Models 1 and 2a. The ratios are 
plotted for two mass ranges of the massive mode: 40-100 M Q (Model 1) and 140-260 M© (Model 
2a). As in the previous figures, the solid line (dashed line) stands for the ISM (IGM). Horizontal 
thin dashed lines indicate the measured abundances in the following four very iron-poor halo stars: 
CS 22949-037 (Depagne et al. 2002), HE 0107-2240 (Bessell et al. 2004), HE 1327-2326 (Frebel 
et al. 2005) and G 77-61 (Plez & Cohen 2005). 
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of increasing U2 by a factor of 5 leads to decrease in Si production by a factor of about 30, with 
little change in C. While it is tempting to argue for a model such as 2a on the basis of [Si/C] (Qian 
Sz Wasserburg 2005), given the large uncertainty inherent in the UV model, it is not clear which 
model (early star formation and chemical evolution or the UV background model) is being tested 
by these observations. The total IGM element abundances inferred from ionized values remain 
uncertain. Nevertheless, Model 1 currently predicts a value of [Si/C] which is too small and may 
signal the need for the presence of some component of more massive stars such as in Model 2a. 



6. Conclusions 

We have calculated the cosmic star formation history corresponding to different minimum 
masses for the initial halo structures, spanning 10 6 to 10 11 M . We include realistic gas outflows 
from the structures, powered by the kinetic energy of supernovae and which take into account 
the increase of the escape velocity due to the growth of structure. We then deduce the baryon 
content and the chemical composition of the structures and of the intergalactic medium (IGM). We 
show that with a minimum halo mass of 10 7 -10 8 M and a moderate outflow efficiency , we are 
able to reproduce both the fraction of baryons in the structures at the present time and the early 
chemical enrichment of the IGM. The intensity of the formation rate of "normal" stars is also well 
constrained by the observations: it has to be dominated by star formation in elliptical galaxies, 
except perhaps at very low redshift. The fraction of baryons in stars is also predicted as are also 
the type la and II supernova event rates. The comparison with SN observations in the redshift 
range z = — 2 allows us to set strong constraints on the time delay of type la supernovae (a total 
delay of ~4 Gyr is required to fit the data), on the lower end of the mass range of the progenitors 
(2-8 Mq) and on the fraction of white dwarfs that reproduce the type la supernova (about 1 
per cent). The type II supernova rate is also well fitted in the same range of redshifts in our 
models and it is directly correlated to the cosmic SFR. We incorporate an improved treatment of 
structure formation compared to our previous work (Daigne et al. 2004) that leads to new insights 
into the initial mass function (IMF) of the population III stars at high redshift. We compare three 
possible mass ranges: 40-100 M Q (normal supernovae), 140-260 M (pair-instability supernovae) 
and 270-500 M . 

We have demonstrated that the fraction in the initial starburst of zero metallicity stars below 
270 M must be limited in order to avoid premature overenrichment of the IGM. Specifically, we 
predict that about 10 - 20 % of the metals present in the IGM at z = have been produced by 
population III stars at very high z. The remaining 80 - 90 % are ejected later by galaxies forming 
normal stars, with a maximum efficiency of the outflow occurring at a redshift of about 5. In full 
agreement with Daigne et al. (2004), because of the chemical constraints, including both for the 
IGM and very metal-poor halo stars, 10~ 3 of the baryons must lie in the first massive stars in order 
to produce enough ionizing photons to allow early reionization of the IGM by z ~ 15. The case of 
the very massive mass range (270-500 M ) is highly efficient regarding the ionizing flux but cannot 
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reproduce alone the global chemical evolution. A massive component of stars with masses in the 40 
- 100 Mq account for the observed ISM and IGM abundances with the possible exception of IGM 
Si. While more massive stars in the range 140 - 260 M produce ample amounts of Si, they can not 
produce the C,0/Fe ratios observed in extremely iron-poor stars, nor are they able to efficiently 
reionize the IGM at high redshift. If future observations bear out a high Si/C ratio in the IGM, 
we would be led to consider a hybrid model of massive stars. 

In summary, we have demonstrated the sensitivity of chemical evolution to different constraints, 
including the minihalo mass range, the role of outflows and the redshift of structure formation, 
correlated with different ranges of massive stars. We have evaluated the rates as a function of 
epoch of both SNII and SNIa. We conclude that only 10 to 20 % of the metals in the IGM are 
produced by Population III. Addition of a massive star component to Population III with masses 
in the conventional stellar mass range provides an effective means of simultaneously accounting for 
both early reionization and the chemical evolution both of the oldest extreme metal-poor stars and 
of the IGM. 
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